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Spatial gradients in refractive index have been used extensively in acoustic metamaterial applica¬ 
tions to control wave propagation through phase delay. This study reports the design and experi¬ 
mental realization of an acoustic gradient index lens using a sonic crystal lattice that is impedance 
matched to water over a broad bandwidth. In contrast to previous designs, the underlying lattice 
features refractive indices that are lower than the water background, which facilitates propagation 
control based on a phase advance as opposed to a delay. The index gradient is achieved by varying 
the filling fraction of hollow, air-filled aluminum tubes that individually exhibit a higher sound speed 
than water and matched impedance. Acoustic focusing is observed over a broad bandwidth of fre¬ 
quencies in the homogenization limit of the lattice, with intensity magnifications in excess of 7 dB. 

An anisotropic lattice design facilitates a flat-faceted geometry with low backscattering at 18 dB 
below the incident sound pressure level. Three dimensional Rayleigh-Sommerfeld integration that 
accounts for the anisotropic refraction is used to accurately predict the experimentally measured 
focal patterns. 

PACS numbers: 43.20.+g,43.30+m,43.40.+s 


Metamaterial lattices composed of sub-wavelength scattering components have been used increasingly in recent 
years to control the propagation of both electromagnetic and acoustic waves through two- and three-dimensions. One 
primary goal of acoustic metamaterial design has been to achieve effective fluid, or “metafluid,” material properties 
that minimize shear coupling and propagation. The coupling between a fluid-born acoustic wave and a composite 
elastic material at off-normal incidence results in mode-mixing of dilatation and shear modes that complicates the 
interaction. For underwater applications requiring strong acoustic coupling, compliant materials such as rubbers have 
been a standard material used for coupling and encapsulation due to their relatively close impedance match with water 
and low shear moduli. Traditional compliant materials often have larger densities and/or lower dilatation moduli than 
water, making them useful for applications requiring low relative sound speed. For example, a sonic crystal lattice of 
rubber cylinders with a gradient filling fraction was recently used to achieve an underwater omnidirectional absorption 
coating [l|, which requires a radially decreasing sound speed. Hi 

Metafiuids with complementary material properties compared to rubbers, e.g. high relative sound speeds with 
matched impedance, have been more difficult to achieve in underwater applications because they require high stiffness 
but low relative density. Sonic crystal lattices feature metafluid functionality in the homogenization limit of the lattice, 
including shear decoupling and broadband performance, and hence they have the potential to expand the range of 
realizable metafluid material properties to include high sound speeds. Thin-walled, hollow elastic shells have been 
recently proposed as high-sound speed sonic crystal components in water. Q By carefully tuning the wall thickness 
of the shell, an impedance-matched condition can be obtained assuming the shell material has higher impedance than 
the background fluid. Additionally, by using shell elements which are individually impedance matched, it has been 
shown that any phononic costal configuration of the elements will also be impedance matched in the homogenization 
limit of the crystal lattice. Q Given that local variations in the element filling fraction create a variation in refractive 
index, acoustically transparent devices with tunable wave-guiding capability should be achievable using hollow shell 
elements. The concept of hollow shell elements has also been extended to include structural components inside the 
shells to improve control over the material properties. [5| 

Here we report the design and experimental demonstration of an acoustically transparent gradient index (GRIN) lens 
that focuses sound by advancing the phase of waves propagating through an aqueous background. GRIN geometries 
have been utilized for a broad range of metamaterial applications including scattering reduction W , wave focusing jsi- 
0 and bending CMiliSlil. In contrast to previous lens designs, we present a lens composed of impedance-matched. 
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hollow shell elements with sound speeds higher than water. The higher sound speed enables the phase of propagating 
waves to be advanced beyond what is possible in the background propagation medium, which is important to a number 
of metamaterial applications (for a review, see Refs.[l3,[^. The lens is intended as a proof of concept to demonstrate 
that both transparency and broadband wave-guiding functionality can be achieved simultaneously using a lattice of 
non-compliant, low refractive index scattering elements. 

The GRIN lens is constructed using cylindrical, air-filled aluminum tubes arranged in a lattice that function as 
a broadband, uniform effective medium at wavelengths A larger than the lattice spacing a (A > 4a). With a wall 
thickness of l/20th the diameter, each individual cylindrical tube scatters sound as an impedance-matched effective 
fluid cylinder Q resulting in negligible backscattering over the lattice as a whole. An anisotropic lattice spacing is 
used to vary the filling fraction similar to what was proposed by Lin, Tittmann, and Huang [1^. The anisotropic 
design simplifies the lens construction by using a single cylindrical tube geometry and features a flat lens facet. While 
the objective of many metafluid applications is to achieve anisotropy in the wave speed, the transparency of the lattice 
results in an approximately isotropic scattering configuration despite the underlying anisotropy of the design. Our 
design goal is to tune the anisotropy to produce a convex refraction gradient that focuses sound while simultaneously 
minimizing lens aberration. 

Figure [TJa) shows the rectangular unit cell of the cylinder lattice. A variable refractive index can be achieved by 
lengthening one of the lattice constants ay compared to the lattice constant ax in the orthogonal direction. Phase 
speeds Cxx and Cyy in the primary Cartesian directions are shown in Fig. [Hb) as a function of the anisotropy ratio 
ay!ax for a fixed aluminum tube outer radius R = O.da^^. The double index on the phase speeds delineates components 
in an anisotropic tensor. Phase speeds c with an over-bar are normalized to the water background, which is assumed 
to have density = 1000 kg/m^ and sound speed = 1480 m/s. The phase speeds are derived from the longitudinal 
dispersion bands of the acoustic band structure calculated for each anisotropy ratio. The anisotropy in phase speed 
can also be derived by considering multiple scattering effects in the lattice. 

Examples of the acoustic band structure are shown in Fig. [Hc,d). Blue and red lines, with slopes corresponding 
to the phase speeds in Fig. Hb), demonstrate that the longitudinal bands have linear dispersion up to an upper 
cutoff frequency ujaxl^iici) ^ 0.3. The cutoff frequency is constrained by the wrapping of the longitudinal band at 
the Brillioun zone boundary in the TY direction for the highest anisotropy ratio considered. The longitudinal band 
in the FM direction in Fig. [TJc) has a slowness that deviates slightly from circularity due to the underlying cubic 
symmetry. In addition to the longitudinal bands, narrow resonance bands arising from the core-shell architecture of 
the unit cell are predicted at various frequencies. [5| While we do not observe evidence of coupling to these bands in 
our measurements, the resonant bands should not, in general, be overlooked depending on the lattice geometry. 

Our lens design consists of a constant lattice spacing in the x-direction, which produces a flat lens facet, while a 
prescribed lattice anisotropy in the ^-direction produces the convex GRIN profile [see Fig.[2][b)]. The primary impact 
of increasing the anisotropy ratio is to decrease the effective sound speed irrespective of propagation direction. As 
is evident in Fig. [Hb), the difference between Cxx and Cyy is small even at ayjax = 2, resulting in approximately 
isotropic acoustic transport over the entire Brillioun zone. Therefore, an approximate, isotropic sound speed Cavg = 
^xx + ^yy) /2 was uscd to initially guide the design of the GRIN geometry. A number of GRIN profiles that result in 
convex leasing are commonly used in the literature; our design combines two common refraction profiles in order to 
minimize lens aberration in the “optical” (ray) transport limit. 
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where h is the lens height in the gradient direction, cq = Cavg(0) and Ch = Cavg{hl2) are the extremal sound speeds 
obtained from the chosen range of ayjax, and /3 is a parameter that mixes the two GRIN profiles. 

An iterative ray tracing routine based on the eikonal approximation [22| was used to determine the optimal value 
of p that minimizes the lens aberration. The isotropic refraction gradient in Eq. ([T]) prescribes Cavg{y) and hence 
the lattice geometry ayjax for a given value of /3; however, to optimize jS an eikonal approximation that accounts 
for the anisotropic sound speeds in the lattice was employed to improve the accuracy of the physical geometry. The 
differential equation governing the eikonal function ^(x) in the presence of anisotropy can be derived for acoustic 
metafiuids in a similar manner to the electromagnetic case, 
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where ^[x,y] = d^/d[x, y]. The off-diagonal sound speed Cxy can be calculated from the band structure or by considering 
multiple scattering effects. 0 The off-diagonal term was found to be negligible compared to the on-diagonal terms 
for each anisotropy ratio in our lattice. 

Figure [2](a) shows the ray paths derived using Eq. (|4]) for a GRIN lens insonified by a plane wave after optimizing 
the parameter /3. The lens has thickness lOa^^, height h = 27.la^^, and lattice spacing ax = 15.6 mm; these parameters 
were also used to manufacture the physical lens. The choice of /3 = 0.65 produces a focal point at x 1.2 m from the 
center of the lens with minimal aberration of the ray trajectories. The /^-optimization routine starts with the average 
sound speed profile Cavg{y) as input and converts to the anisotropic phase speeds from Fig. [IJb). The optimization 
routine also ensures that the gradient in the lattice spacing ay{y) is properly discretized to be consistent with the 
total height h of the optimized GRIN profile. 

A schematic of the upper half of the /^-optimized lens design is shown in Fig. EJb). The full lens is symmetric 
about the ^-axis. A photograph of the assembled lens is shown in Fig. [2fc). The lens is constructed of identical 1 m 
long hollow aluminum cylinders with air inside. Each cylindrical tube has a diameter of 12.5 mm and nominal wall 
thickness of 0.625 mm. Rubber end caps are adhered to both ends of each cylinder to prevent water infiltration. 
Additionally, ^2.5 cm of foam was inserted into the tube ends to damp the excitation of axial modes. The tubes are 
mounted on plexiglass end-plates into which holes with the designed lattice pattern are drilled. Two additional rows 
of cylinders with spacing ay = ax are placed at the top and bottom of the lens to help mitigate the sharp transition 
where the edge of the lens {cavg = 1.36) meets the open water {cavg = !)• 

The lens was submerged at the center of a 6 x 6 x 4 m^ water tank. Acoustic waves were produced by a 0.1 m 
diameter spherical source located in the plane that bisects the lens midpoint in the axial direction ( 2 ;-axis). In order 
to demonstrate the lens directionality, two in-plane source locations were considered at {x^y) = (—1.6,0) m and 
{x^y) = (—1.6,0.43) m corresponding to incident angles of 0° and 15° respectively from the lens central axis (x-axis). 
Wave propagation was measured using hydrophones at a sampling rate of 0.8 MHz. Hydrophones were mounted on a 
three-axis translation positioning system to record measurements of the sound pressure levels, P(x, y) and Po(^, y)i in 
the presence and absence of the lens respectively. Measurements were taken at 10 mm increments. The sound pressure 
level was measured at individual frequencies by averaging over a 10-cycle pulse in the transmission region (as indicated 
in Figs. [3] and S]) and a 20-cycle pulse in the refiection region. The longer pulse cycle used in the refiection region 
ensured overlap between incident and refiected waves so that the measured total field could be properly compared 
with simulations. The length of the incident pulses were short enough to isolate refiections from the tank walls. 

Measurements were performed in the x^-plane perpendicular to the cylindrical axis of the lens over a broad range 
of frequencies ujaxl^iici) < 0.3. Examples of the measured and numerically modeled pressure intensity (|P|/|Po|)^ in 
the vicinity of the lens are shown in Figs. [3] and |4] for frequencies ujaxl^iici) = 0.077, 0.138, and 0.256. In both the 
measured and predicted results, the pressure intensity is normalized (at each position) to that of the source intensity 
in the absence of the lens. Given the symmetry of the lens about the ^-axis, measurements were obtained over 
the bottom half of the scattering plane with a small overlap into the upper half-plane to detail the on-axis forward 
scattering pattern. The regions scanned in the measurement are also outlined in the upper panels of Figs. [3] and |4] for 
ease of comparison with the numerical predictions. Note that although the measurements at ujaxl^'KCij = 0.077 and 
0.256 lie within two of the resonance bands identified in Fig. [U no obvious additional resonant features are observed 
in the intensity maps that can be attributed to these resonance bands. 

The upper panels (a-c) of Figs. [3] and [4] show a prediction of the transmitted and refiected pressure intensity 
calculated using two-dimensional multiple scattering theory (2D-MST). 0,0 The MST accounts for elastic scattering 
in the cylindrical tubes and assumes insonification by a cylindrical monopole source located at positions that match the 
experiment. The experimentally measured acoustic intensities in panels (d-f) show good qualitative agreement with 
the MST-predicted results. The refiected signal, R = 201og]^o[(^ “ Po)/Po)]i was estimated relative to the incident 
amplitude Pq using the sound pressure level P measured in front of the incident face of the lens (regions to the left of 
the lens in Figs. [3] and 0]). The measured refiection was at or below 13% of the incident amplitude {R < —18 dB) over 
the range of operational frequencies, which demonstrates significant transparency and is commensurate with recent 
Fresnel lens designs in air. [^, The measured refiection includes additional constructive interference compared to 
the MST due to diffraction from the finite aperture size in the axial direction. 

In the forward direction a focusing peak is observed that strengthens in intensity and moves out away from the lens 
with increasing frequency. The frequency-dependence of the focusing peak is expected. At very low frequencies the 
aperture is diffraction-limited; the wavelength is too long to adequately resolve the index gradient and there is also 
significant diffraction around the aperture. As the frequency increases, the transmission approaches the geometric 
limit of ray-acoustics where a focusing peak would be located beyond the focal point predicted by the ray tracing in 
Fig. [2ja). At the lowest example frequency, ujaxl^nci) = 0.077 [panels (a,d)], the lens is close to the diffraction limit 
with height-to-wavelength ratio h/X < 2. As the frequency is further reduced the focusing peak intensity becomes 
significantly suppressed. Therefore, an operational bandwidth can be identified ranging between uaxl2'KC\y ~ 0.05- 
0.30, where the lower cutoff is constrained by the diffraction limit and the upper cutoff by the homogenization limit 


4 


of the lattice. 

It is apparent in Figs. [3] and 0] that the 2D-MST does not accurately predict the precise location of the intensity 
maxima nor the small-scale, near-held interference features in the forward direction. The discrepancy is due to the 
assumption of inhnite cylinders in the 2D theory compared to the hnite size of the lens aperture in the experiment. 
The diffractive corrections arising from the hnite size in the axial direction can be predicted using a three-dimensional 
(3D) Rayleigh-Sommerfeld integral. The pressure at a position (x, y, z) in front of an acoustic aperture located at 
X = 0 can be calculated by integrating over the aperture, [25| 

r r ikT" 

P{x,y,z) = ^ jj P(y'^z'f—dy'dz' (5) 

A 

P{y',z') = (6) 

= A + (y - y'f + {z- z'f (7) 

where and (j){y'^z') are the real-valued acoustic amplitude and phase of the complex pressure P{y'^z') on 

the aperture facet, and A and k are the wavelength and wave number respectively. Here the integration aperture is 
simplihed by the hat facet of our lens design, with the integration carried out over a simple rectangular area at the 
forward face of the lens. 

The pressure on the aperture facet P{y'^z') can be calculated using anisotropic ray-tracing under the eikonal 
approximation of Eq. ( 0 ). A schematic in Figure [5fa) depicts the ray trajectories of wave-fronts that emerge from the 
spherical source and traverse the forward facet of the lens. The calculation iterates over a discrete set of rays that 
enter the rear facet (i.e. the incident face) at equally spaced intervals. The phase function (j){y'^z') on the forward 
facet was determined by advancing the phase of each ray along the path lengths defined by ^(x) assuming a constant 
phase at the source. The effective path length is altered by the spatially-dependent local sound speed in the lens. 
The phase difference A0 = (j){y'^ z') — 0(0, 0) relative to the phase at the lens central axis is shown in Fig. Efb), where 
positive or negative values represent a phase advance or delay respectively. The calculation is performed for a source 
position on the lens central axis. Given the symmetry of the lens only one quadrant is shown. Curve fits in Fig. [5fd) 
indicate that A0 is approximately quadratic in both the y- and 2 ;-directions. 

The amplitude function 0(^0 z') is inversely related to the square root of the local areal density of the rays. This was 
approximated by calculating the average nearest-neighbor separation of rays intersecting the forward facet; a relative 
amplitude fj{y'^z') = fj{y', z')/fj{0,0) can then be estimated from the inverse of this nearest-neighbor separation. 
Figure [5](c) shows the variation in {^{y'^z') after normalizing to the predicted amplitude in the absence of the lens. 
It is clear that the lens refraction has minimal impact on the amplitude (< 3%) and that the diffractive correction is 
primarily influenced by the change in phase over the aperture. 

Following a similar procedure as was used by Gao et. a/., [26[ Eq. (j5j) is evaluated numerically within the 2D plane of 
the measurement using quadratic approximations to 0(^', z') and 0(^0 z'). Intensity maps produced by the Rayleigh- 
Sommerfeld integration are shown in Fig. [3l^g-i) and demonstrate significantly improved quantitative agreement with 
the experiment. As was the case for the measurement and MST modeling, the Rayleigh-Sommerfeld-based pressure 
intensities are normalized to the pressure intensity of the spherical source in the absence of the lens. The 3D Rayleigh- 
Sommerfeld integration more accurately predicts both the location of the focal positions and their magnitude. The 
measured magnification in sound pressure level ranges between 4.5-7.3 dB over the operational bandwidth. The 
2D-MST predicts magnifications that are as much as 40% lower than the measured values, whereas the magnifications 
predicted by the Rayleigh-Sommer feld integration agree to within 10%. We emphasize that the 2D lattice only focuses 
along one axis (^-axis) resulting in lower magnifications than would be produced by axisymmetric 3D designs. Hilil 
The magnification can be increased by extending our 2D design to three dimensions using tubes bent into a toroid 
configuration. mui 

We now place our results in context with other metafluid design concepts. Traditional underwater and ultrasound 
applications have achieved an impedance-match by utilizing rubbers that have sound speeds less than water. While 
rubbers have been used as components in transparent metafiuid lattices, [I| there have been other recent advances in 
acoustic impedance-matched designs. For example, negative-index complementary metamaterials have been proposed 
to significantly enhance the transparency through aberrating materials. A class of metafluids based on the concept 
of space-coiling have also featured prominently in the literature. Space-coiled Fresnel lens desi gns have been 

reported that show significant focusing while utilizing a thin aperture compared to the wavelength. [^, Although 
the strength of the space-coiling design is the ability to significantly modify phase within a confined space, there is 
a drawback that it can only delay the phase over a lengthened path; similar to rubbers in water, the effective sound 
speed of these devices is less than the background fluid. The space-coiling design has not yet been demonstrated in 
water, where a higher viscosity and lack of rigid boundary conditions may impose additional challenges. 

Our design extends the reach of aqueous transparent metafluids to include the option of phase advanee and decreased 
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effective path length. While not necessarily required for traditional lensing, there are metamaterial applications 
that require higher sound speeds compared to the propagation medium, |17l - [^. 1^ the most prominent of which is 
scattering reduction. The requirement of higher relative sound speeds features prominently in scattering reduction 
designs based on both coordinate transformation j32l-[^ and scattering cancellation. [sH, [36| The realization of a high- 
sound speed, transparent GRIN lens represents an important proof-of-concept: it demonstrates that phase-advance 
metafluids can be constructed with significant tunability in both sound speed and impedance, including the option of 
matched impedance. Given that the concept of core-shell lattice components has now been extended to also include 
core modification, [5| a wide range of GRIN and anisotropic designs should be achievable in water based on metafiuid 
lattices. 

[Work is supported by the Office of Naval Research.] 
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FIG. 1: (a) Schematic unit cell of the anisotropic lattice with ay = 2ax. (b) Phase speeds in the homogenization limit of the 
lattice plotted as a function of anisotropy ratio ayjax- Lines are hts to the data based on Ref. 0. (c) Acoustic band structure 
calculated using the Finite Element Method (FEM) for an isotropic lattice with tube outer radius R = OAax. Blue lines plot 
the acoustic sound speed in the FX direction, (d) Acoustic band structure calculated for an anisotropic lattice with tube outer 
radius R = OAax. Blue and red lines plot the acoustic sound speed in the FX and FY directions respectively. Insets show the 
directions of the bands in /c-space. 
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FIG. 2: (a) Ray paths through the lens (blue lines), which originate from a plane wave at normal incidence, are calculated 
using an eikonal approximation in the xy-plane; the red box indicates the location of the lens. Only rays incident on the upper 
half-plane of the lens are shown, (b) Schematic showing the locations of the lattice sites in the upper half-plane of the lens; 
the resulting sound speed profile is plotted to the right. The lattice sites are symmetric about y = 0. (c) Photograph of the 
assembled lens composed of 1 m long, hollow aluminum cylinders arranged in the desired lattice pattern. 



FIG. 3: Sound pressure intensity in the vicinity of the lens at uaxl^Tvcb = 0.077 (7.27 kHz), 0.138 (13.05 kHz), and 0.256 
(24.24 kHz) for a source location on the x-axis. (a-c) Pressure intensity calculated using 2D-MST. (d-f) Measured pressure 
intensity, (g-i) Pressure intensity calculated using a 3D Rayleigh-Sommerfeld approximation. Grey boxes show the position of 
the lens. Red outlines indicate the experimentally mapped areas. The reduced observational range of the reflected signal at 
24.24 kHz is due to a limited overlap between the incident and reflected signals in the time domain at high frequency. 
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FIG. 4: Sound pressure intensity in the vicinity of the lens at uaxl^Tvcb = 0.077 (7.27 kHz), 0.138 (13.05 kHz), and 0.256 
(24.24 kHz) for a source positioned at a 15° angle with respect to the a:-axis. (a-c) Pressure intensity calculated using 2D-MST. 
(d-f) Measured pressure intensity. Grey boxes show the position of the lens. Red outlines indicate the experimentally mapped 
areas. 
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FIG. 5: (a) Schematic showing the experimental setup; blue lines indicate ray paths emerging from the spherical source and 
Rayleigh-Sommerfeld integration vectors from the forward facet of the lens to hydrophone locations in the measurement plane, 
(b) Change in phase A0(y', z') (degrees) over the forward facet (due to symmetry only the 1st quadrant is shown). Blue regions 
represent a phase delay, red regions represent a phase advance, (c) Normalized amplitude '0(y^ z') over the forward facet, where 
normalization is with respect to the amplitude of the spherical wave in the absence of the lens, (d) Red and blue lines show 
the phase change A0 plotted as a function of position along the y'- and z'-axes respectively. Data points are calculated by 
advancing the phase along ray paths based on an eikonal approximation, lines are quadratic fits to the data. 



















